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ABSTRACT 

We present the results of a search for gravitationally-lensed giant arcs con- 
ducted on a sample of 825 SDSS galaxy clusters. Both a visual inspection of 
the images and an automated search were performed and no arcs were found. 
This result is used to set an upper limit on the arc probability per cluster. We 
present selection functions for our survey, in the form of arc detection efficiency 
curves plotted as functions of arc parameters, both for the visual inspection and 
the automated search. The selection function is such that we are sensitive only 
to long, high surface brightness arcs with g-hand surface brightness /ig < 24.8 
and length-to- width ratio l/w > 10. Our upper limits on the arc probability 
are compatible with previous arc searches. Lastly, we report on a serendipitous 
discovery of a giant arc in the SDSS data, known inside the SDSS Collaboration 
as Hall's arc. 

Subject headings: clusters, arcs, SDSS 



1. Introduction 



Clusters of galaxies are the largest gravitationally bound structures in the universe. It 
has long been recognize d that they constitute very us efeul probes for cosmology (see, e.g.. 



Press fc Schechterl Il974j : iMohr et al.l l200ll : IVoitI l2005l ). and therefore understanding them 
in detail has been a major field of research. Notable progress in the theoretical under- 
standing of clusters has been achieved in recent years through N-body and hydrodynamical 
simulations. At present, different simulation schemes converge in reproducing several broad 
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properties of clusters (iFrenk et al.lll999l : iHeitman et al.l l2005l ) . Furthermore, several prop- 
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The strong gravitational lensing of background galaxies produced by massive clusters is 
a probe of cluster structure. This effect has been observed in man y previous studies, w i th the 
mos t extreme case being the giant gravitational arcs first seen by iLynds fc PetrosianI (119891 ) 
and ISoucail et al.l (119871 ). Almost a decade later, spectacular image s of these arcs were 



provided by the excellent resolution of the Hubble Space Telescope (I Couch fc EUid Il995 



Kneib et al. 


1996; 


Small et al. 


1996) 



There are n ow about a hundr ed clusters with giant 



arcs imaged by HST (for a recent compilation, see ISand et al.ll2005l ). Several arc searches 
have been carried out with ground based o 
e.g 



Luppino et aP Il999l: Icvprianol l20oi 



jservations, both for X- ray selected clusters (see. 



Campusano et al.l 



identified clusters ( Zaritsky fc Gonzalej 20031 : Gladders et al. 



200 



^ as well as for optically 



20031). 



The comparison of the abundance of gravit ational arcs with theoretical predictions has 
been proposed as a tool to c onstrain cosmology (iGrossman fc Narayanlll988l : iMiralda-Escude 
19931 : IWu fc Hammerlll993l ). The early results indicated that the number of strongly -lensed 
arcs greatly exceeded that expected from ACDM simulations (IBartelmann et al.lll998l ). More 
detailed studies of the theoretical uncertainties and systematic effects involved in this analysis 
seem to indicate that the obser vations do, in fact, agree with the theoretical expectations 
from the standard cos mology ( Hennawi et al. 2005 : Dalai et al. 2004 : ^Wamb sganss et al. 



2004J : lOguri et al.ll2003l ). although uncertainties may still remain (ILi et al.ii2005l) . 



The main goal of the work presented here is not in the extraction of cosmological infor- 
mation from the statistics of gravitational arcs. Instead it is both (a) to locate high surface 
brightness giant arcs that would prove useful in follow-up imaging and spectroscopy pro- 
grams, and more directly (b) to understand the probability for a galaxy cluster to produce 
a giant arc as a function of the cluster mass and redshift. We believe this is an important 
aspect to be considered before cosmology can be extracted from arc statistics and that it 
will provide a significant improvement in our understanding of the mass distribution inside 
galaxy clusters. 

A defini ng characteristic of pre vious arc searches is that they were done in the most mas- 
sive clusters. iLuppino et al.l (119991 ) searched for arcs in 38 X-ray selected clusters (redshifts 
0.15 < 2; < 0.82) with X-ray luminosity >2 x 10^"^ ergs/sec and found arcs in a high frac- 
tion of clusters with > 1 0"^^ ergs/sec, but no arcs in clusters with < 4 x 10'^^ ergs/sec. 
Zaritsky fc Gonzalea (120031 ) searched for arcs in 44 optically selected c lusters (0.5 < z < 0.8) 
from a list ranked by surface brightness, finding arcs in two clusters. iGladders et al.l (120031 ) 
searched ~ 900 optically selected clusters (0.3 < z < lA) ranked by a galaxy overdensity 
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parameter, finding arcs in 8 clusters, with none ai z < 0.64. ISand et al.l (120051 ) used a very 
different technique, searching for arcs in the heterogeneous sample of clusters that have been 
observed with HST; out of 128 clusters (0.1 < z < 0.78) 45 had tangential arcs, although it 
is worth noting that many of the clusters were targeted by HST precisely because of known 



arcs. 



It is possible to compute which of the arcs found in high quality images would be found 
in images of lower quality, and we will perform such a comparison. Th e results so far suggest 
looking for arcs in clusters with high X-ray luminosity or high redshift. iLuppino et al.l (119991 ) 
found a higher frequency of arcs in high clusters, and since the bremstrahlung emission 
in clusters is proportiona l to n^, this resu l t sugg ests that arcs are found preferentially in the 
highest density clusters. iGladders et al. J2003h found a higher frequency of arcs in higher 
redshift clusters, which suggests either an evolution of cluster structure in a way not fully 
understood, or perhaps that the sources lensed are predominately at higher redshifts, where 
the lensing cross se ction would be much larger for z ~ 0.7 clusters than for those at 2; ~ 0.3 
(IHo fc Whitd I2OO5I ). Observationally, however, it is clear that the dominant variables for 
successful arc detection are limiting surface brightness and, in particular, seeing. 



The Sloan Digital Sky Survey (lYork et al.ll2000l ) to date has produced multiband ugriz 
imaging of 8000 deg^ of the norther n sky, along with an associated catalog of 1,048,496 
spectra, of which 674,749 are galaxies (lAdelman-McCarthy et al.ll2006l ). It has proven useful 
as a data set for gravitationa l len s searches. An ex ample is the multiply-lensed quasar 
searches of iPindor et al.l (120031 ) and lOguri et al.l (120061 ) . which use the imaging data to select 
quasar- colored objects that are larger than a PSF, or groupings of quasar-colored objects 
in a small area. A sec ond example of lens searches in the SDSS data uses the large SDSS 
spectroscopic data set: iBolton et al.l (120061 ) searched red galaxy spectra for indications of a 
second, higher redshift nebular emission spectrum. This locates strongly-lensed background 
galaxies, which are revealed as arcs after applying galaxy subtraction techniques on HST ACS 
images. Our program is a search for strongly-lensed background galaxies behind clusters; 
that is, a search for giant arcs. Our ability to locate gravitational arcs in the SDSS images 
is demonstrated by recovering arcs observed by other surveys and by the discovery of new 
gravitational arcs (see Fig. [T|, Appendix \^ and lAUam et al.l 120061). A simil a r arc search in 
a sample of 240 rich SDSS clusters has been carried out by iHennawi et al.l (120061 ) and has 
found a significant number of new giant arc systems, but relies on deeper imaging data to 
discover the arcs, whereas we use the original SDSS imaging data itself. 

Our analysis consists of a search for giant gravitational arcs in a sample of SDSS galaxy 
clusters sorted by a richness estimator. Our imaging data does not have the depth or seeing 
available to previous searches, but it does have the twin features of homogenity and large 
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sky coverage. This results in our ability to greatly increase the number of clusters searched. 
The search reported here totals 825 clusters, comparable to the largest previous search. We 
intend, in a future paper, to decrease the mass threshold down to the poor group range in 
order to search tens or hundreds of thousands of locations. Furthermore, the redshift range 
of the cluster catalog we have searched is predominantly 0.1 < 2; < 0.3, a regime poorly 
covered by previous searches. 
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Fig. 1.— On the left, the SDSS image of the RCS141 9.2+5361 cluster. Th e arc was dis- 
covered by visual inspection of the RCS cluster sample (iGladders et al.ll2003l ) and is clearly- 
visible in this color image composed of the SDSS r, and i data. On the right, Hall's 
arc (see Appendix), serendipitously discovered in the SDSS data. It too is clearly visible in 
the composi te qri image. The se observations, together with the recent discovery of a new 
arc system (lAUam et al.ll2006l ). demonstrate the ability of the SDSS to detect high surface 
brightness giant arcs. 



- 6 - 



The Data 



Adelman-McCarthv et al. 


2006) 


camera ( 


Gunn et al. 


1998 


) on 



The SDSS uses a wide-fiel d drift-scanninK m osaic CCD 



Point Observatory, New Mexico, to image the sky in the five photometric bandpasses ugriz 



( jFukugita et al. I Il996l : IStoughton et al.ll2002l ). The SPSS photomet ric data are processed 



using pi pelines that address astrometric calibrations ( Pier et aD 2003h. photornetric dat 



duction (ILupton et al 



20011). and photom etric calibrations (ITucker et al.l l2006l : ISmith et al. 



ca re- 



2OO2I : iHogg et al.ll200ll : llvezic et al.ll2004f ). Note that the SDSS is a survey that trades rela- 
tively shallow exposures - 55 seconds on a 2.5m telescope, resulting in point sources detected 
at 95% confidence level at g = 22.2 and z = 21.3 (AB) - for a very wide area: 8000 deg^ or 
20% of the total sky. 



2.1. The Cluster Sample 



The cluster catalog was constructed using the maxBCG algorithm. This is a red- 
sequence locating technique that has shown completeness levels of 90% at 0.1 < 2; < 0.3 
and down to A^^^ = 10, where Ng is the number of galaxies on the E/SO red sequence brighter 
than 0.4L* and within lh~^ Mpc of the brightest cluster galaxy. The correlation of the 
richness Ng with mass can be calibrated with independent mass estimators, such as velocity 
dispersions or weak lensing shear pro fil es. A detailed descrip tion of the maxBCG technique 
is discussed in iKoester et al.l ( l2006al ). iKoester et all (l2006bl ) used the SDSS DR5 data to 
construct a cluster catalog containing 12,875 rich clusters, plus two orders of magnitude 
more clusters of lower masses. The catalog we used in this work was created using an earlier 
version of the maxBCG algorithm (see, e.g.. lHansen et al.ll2005l : ISheldon et al.ll200ll ). 



One of the defining features of red sequence methods is that projection effects are present 
only at a minimal level: if the clusters differ by Az < 0.05 then they can be confused as a 
single object. The resulting purity (lack of false positives) is ~ 95%, and the completeness 
(lack of false negat ives) is ~ 85% above 1 x 10^"^ Mq, as determined by running this algorithm 
on mock catalogs (IKoester et alll2006bl ). 



We selected a sample of 825 clusters with Ng > 20, approximately 50% of the Ng > 20 
clusters in our sample. The whole sample was not searched due to a variety of technical 
reasons, but lack of a complete sample is not important for the results of this paper. The 
distribution of (yf-band seeing of the cluster images is shown in Fig. [2] and has a mean of 
1.44"; this is important because the seeing is the single dominating factor in our ability to 
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detect arcs. The cluster redshift distribution is shown in Fig. [31 The richest clusters in our 
sample are presented in Table O and the table listing our full cluster sample is available 
electronically. 
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Fig. 2. — The g-hcuid seeing distribution for the SDSS images of the 825 clusters inspected 
in this paper. The distribution is well fit by a gaussian of mean 1.44" and standard deviation 
0.20". 
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Fig. 3. — The redshift distribution for the 825 clusters inspected for this paper. We have 
good statistics in the redshift range 0.1 < 2; < 0.3. 
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While we do not use the masses of the clusters directly in this work, we do have a rough 
mass cahbration. Stacked cluster spectroscopic velocity dispersions, plus isothermal sphere 
fits to stacked cluster weak lensing shear profi les, both from early work in the SDSS EDR (see. 



e.g 



Bahcall et al.l l2003l : ISheldon et al.ll200ll ). were combined to produce a scaling relation 
between velocity dispersion and Ng given by ~ 85iV°-^ km/s. Modehng the cluster as 
a singular isothermal sphere, the mass enclose d within a radius having an overdensity of 200 
with respect to the critical density is (see, e.g. iBryan fc NormanI 119981 ) 



200 



1 

a J E{z)h 



0) 



where (Tq = 1148 km/s, E{z) = ^/Qa + (^^(1 + z)^ + {1 - Qa - f^m)(l + z)^, h is the 
Hubble constant in units of 100 km/s/Mpc, and Qm and are the energy densities of 
cold/presureless normal matter and the cosmological constant, respectively, in units of the 
critical density. The a^-Ng relation has a significant scatter. Indirect methods, such as (a) 
running the cluster finding algorithm on mock catalogs based on N-body simulations and 
(b) measuring the higher moments of the velocity distributions of galaxies in stacked cluster 
samples, suggest that the scatter is AM2oo/^2oo ~ 40% at high Ng and increases at lower 
Ng. A more sophis ticated calibrat ion u sing weak lensing mea surements of the full maxBCG 



catalog is given in ISheldonI (120061 ) and IJohnston et al.l (120061 ). In future work we expect to 
use the current maxBCG catalog and mass calibrations. 
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Fig. 4. — The mass distribution for our sample of inspected clusters. The masses were 
estimated from the a^-Ng relation mentioned in the text and from Eqn. ([1]). 
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3. Search by Visual Inspection 

Previous arc surveys have generally used visual inspection to locate arc candidates, 
which are then usually followed up with deeper imaging in better seeing. 

The SDSS images of the 825 selected clusters were visually inspected for arcs. For each 
cluster, we used a SDSS coaddition code to build an image covering an angle corresponding 
to a physical scale of 1 Mpc (assuming Hq = 100 km/s, Qm = 0.3, and = 0.7). The 
inspector was then presented with 4 simultaneous images using the ds9 image display pro- 
gram: grayscale images of the g, r, and i bands and a color image that combines the three 
bands. The grayscale images were displayed using the ds9 "zscale" algorithm, while the 
color image was displayed using a fixed surface brightness range (the same range was used 
for all clusters). A single author of this work (T.L.) scanned the 825 images, while three 
other authors inspected images of the 300 highest A^^ clusters in the sample. Candidates 
were recorded for later consideration and for potential follow-up observations (§S])- 



3.1. Visual Inspection Selection Function 

A selection function is a prerequisite for the proper use of arc statistics to constrain 
cosmology. In order to calculate a selection function for our arc search, simulated arcs were 
added to a fraction of the images in our cluster sample. These objects were included before 
the images were inspected for arcs and without the knowledge of the scanner doing the visual 
inspection. Each added arc consisted of a section of a circle centered on the brightest cluster 
galaxy identified by the maxBCG algorithm. We added the arcs using uniform distributions 
in log peak surface brightness and in length-to- width ratio {l/w = 5 — 20) in order to derive 
the arc detection efficiency as a function of these arc parameters. The simulated arcs have 
a flat spectrum in AB magnitudes, similar to typical observed gravitational arcs from other 
surveys (this makes the simulated arcs much bluer than the cluster galaxies). The typical 
length for the simulated arcs was I ~ 10", and they were inserted in the images at 5" — 45" 
from the brightest cluster galaxy (BCG). The simulated arcs have a surface brightness profile 
that is uniform along the tangential direction and is a gaussian along the radial direction, 
with a FWHM = 1.2" (si w), similar to the typical SDSS seeing. The selection functions 
obtained with this technique are presented in Figs. [5] and [61 
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Fig. 5. — Efficency for finding simulated circs clS db function of the length-to-width ratio l/w. 
The panels on the left show the number of simulated arcs included in the images (dashed 
black) and the number of these arcs recovered by the visual inspection (solid red). The panels 
on the right show the efficiency for recovering the simulated arcs. The top panel shows the 
results for arcs with more than 5 counts/pixel in peak surface brightness (the conversion of 
counts per pixel to surface brightness is shown in Fig. [H]) and the bottom panel shows the 
efficiency for the whole sample. 
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Fig. 6. — Efficiency for finding simulated circs cls cL function of gf-band peak surface brightness, 
in units of either counts per pixel (top panels) or AB magnitudes per arcsec^ (bottom panels). 
The panels on the left show the number of simulated arcs included in the images (dashed 
black), the number recovered by the visual inspection (solid red), and the number recovered 
by the automated search (dotted blue). The panels on the right show the efficiency for 
recovering the simulated arcs. 
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For the visual inspection, the detection efficiency reaches 50% at about 24.8 mag/arcsec^ 
in the peak g-band surface brightness and for l/w ~ 10. The primary causes for missed arcs 
are low surface brightness and small length-to-width ratio. However, the recovery efficiency 
does not reach 100% even for high l/w ratio, high surface brightness arcs, due to blending 
problems with bright galaxies and stars on the images. Note that the measurement of our 
selection function is somewhat limited by small-number statistics, because we decided not 
to populate a significant fraction of our images with simulated arcs (we added 65 simulated 
arcs in total, i.e., to about 8% of the 825 inspected clusters). 



4. Search by Algorithm 



While sophisticated arc-finding algorithms have been published previou sly (e.g. , iLenzen et al 



2004 : iHoresh et al.ll2005l : lAlardll2006l : IStarck et al.ll2003l : ISeidel et al.ll2006f ) we are not aware 
of an arc search that has been performed in a fully automated way. Clearly if we are to ex- 
amine 10^ clusters we will need to use an algorithm (as opposed to a human eye inspection). 
We have developed such an algorithm, based on ideas taken from related techniques used in 
experimental particle physics for tracking particles in cloud and bubble chambers. In this 
work we will set a benchmark comparison between a simple algorithm for arc detection and 
a visual inspection. 



4.1. The Algorithm 



Our algorithm is composed of three steps: image preparation, candidate pre-selection, 
and final arc selection. 

Step 1 is image preparation, where we scale the 5',r, i images by the variances of their 
respective sky background s, normalize the variance s to 1.0, average the three images together, 
and then run SExtractor (IBertin fc Arnoutslll996l ) to produce an object image. This step is 
intended to enhance the signal-to-background ratio of the arc-like obje cts. The averaging of 
the variance- normalized images follows the ideas of ISzalay et al.l (119991 ). and the intent is to 
produce an image where the pixel values are a of the hypothesis that the pixel consists 
of only sky. The SExtractor-produced object image has eliminated isolated pixels above 
threshold (due to noise) that do not belong to any object, retaining only pixels belonging 
to objects dectected above threshold. We perform a deblending of the SExtractor objects 
by separating each object into sub-objects which correspond to groups of contiguous pixels 
above threshold, and those sub-objects that are separated from another sub-object by fewer 
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than 2 pixels are merged into one object. 

Step 2 is candidate pre-selection, where the basic shape of the candidates is characterized 
by measuring the major and minor axes of each object, and where objects that show elongated 
shapes are kept for further study. The major axis (ri) is defined as the largest distance 
between any pair of the N pixels belonging to the object, and the minor axis (r2) is defined 
as the maximum distance between pixels perpendicular to the main axis. Elongated objects, 
where ri/r2 > 1.4, are chosen for further analysis. Step 2 simply reduces the number of 
objects to be studied, in order to reduce computation time per image. It is designed to be 
both complete and efficient, keeping 100% of the arc-like objects while providing about a 
factor of 100 in rejection against the rest of the objects detected by SExtractor. 

Step 3 is final arc selection, where the remaining objects have their radius of curvature 
measured. Those radii, plus other parameters measured in step 2, become the input for a 
neural net (NN) trained to select simulated arcs. The output of this NN is used to determine 
if an object is a good arc candidate. The radius of curvature (R) is measured using a least 
squares fit. The selection of good arc candidates then uses the following quantities: 

• Ti/ {27cR): the fraction of the circumference covered by the arc. This is typically > 0.2 
for objects with significant curvature. 

• N/IttR"^): the fraction of the circle area covered by pixels above threshold. This is 
typically < 0.15 for good candidates. 

• /N: the goodness of fit. 

• Ti. the length of the long axis. This ensures significant size for the candidates, > 20 
pixels (1 pixel = 0.396"). 

The 4 quantitites described above are then presented to a NN with 8 nodes in a hidden 
layer and the output is trained (by back propagation) to identify simulated arcs. The NN 
was trained using a separate sample of 100 bright simulated arcs, where the peak surface 
brightness was fixed at 20 counts per pixel in the g band. 

4.2. Algorithm Selection Function 

The efficiency obtained with this algorithm is somewhat lower than that seen for the 
human scan and is shown in Fig. [6l The detection efficiency reaches 40% at 24.7 mag/arcsec^ 
in the g band and for //w ?5i 10; the curve is consistent with the eye scan curve shifted to a 
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brighter surface brightness by ^ 0.4 mags. More development will be needed to bring this 
algorithm to the level of efficiency of the visual inspection. However, this analysis provides 
a quantitative benchmark for the comparison of a human search with a simple algorithm 
search applied to the same cluster sample. 



5. Candidate Follow-up 

The six best candidates, all from the visual inspection, were selected for follow-up. These 
arc candidates were observed using the SPIcam imager on the Apache Point Observatory 
3.5-meter telescope. Images were takin in the SDSS r, and i filters. The observations were 
done the night of April 29-30, 2006, during dark time. Table [1] shows the amount of time 
spent observing each object with the three filters. The images were reduced using standard 
bias subtraction and fiat-fielding techniques. The images were also registered and coadded. 
Typical seeing on the individual images ranged from 1.1" to 1.5". 

The candidates were inspected by the authors and discarded based on judgements that 
they were most likely "junk" or edge-on galaxies. Room for improvement in this system exists 
and we plan on putting into place a set of criteria for systematic follow-up. The criteria useful 
for SDSS data are obvious from an examination of Fig. [TJ l/w, surface brightness, radius of 
curvature, and color. 
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Table 1. Candidate Observations 



Name 


RA 


Dec 


z 


g (min) 


r (min) 


i (min) 


SDSS+139.5+51.7+0.24 


09:17.59.5 


51:42:06 


0.24 


15 


5 


5 


SDSS+184.4+3.6+0.08 


12:17:28.0 


03:36:17 


0.08 


15 


5 


5 


SDSS+117.7+17.7+0.19 


07:50:48.9 


17:40:43 


0.19 


15 


5 


5 


SDSS+175.3+5.8+0.12 


11:41:13.5 


05:48:28 


0.12 


15 


5 


5 


SDSS+205.5+26.4+0.10 


13:41:49.1 


26:22:25 


0.10 


15 


5 


5 


SDSS+139.1+5.9+0.14 


09:16:21.3 


05:53:18 


0.14 


15 


5 


5 



Note. — SPIcam observations. Exposure times are given in minutes. The redshifts 
quoted for the clusters are photometric redshifts from the maxBCG algorithm and 
have a dispersion of cr^ = 0.015. 
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Results of the Searches 



Our survey of 825 SDSS galaxy clusters resulted in no convincing candidates for giant 
gravitational arcs. 

The selection function we have computed shows that we are > 50% efficient at SB < 24.8 
mags/arcsec^ and for l/w > 10. The latter, given our mean seeing of 1.4", corresponds to 
arcs longer than 14". Fig . [7] shows the surface br ightness distribution of the arcs found by 
Lvnds fc PetrosianI (Il989l. in seeing of ~ 1.0"), bv lLuppino et al.l (Il999l. ~ 0.90". their Table 
2). bv lZaritskv fc GonzaB (|20oJ ~ 0.55"), and bv bladders et al.l J2003I . ~ 0.85"). These 
are all the published ground based surveys, and no attempt to homogenize the seeing has 
been made. The surface brightnesses have been converted into the AB system from the 
Vega system. We assume that the arcs have flat spectra so that no change is necessary for 
the different bandpasses. We see from Figure [7] that our surface brightness limit is in fact 
exceeded by the majority of these arcs. The observed surface brightness for arcs that are 
unresolved in width is a function of the seeing obtained, and the data here include a variety 
of seeing. However, all are better than our sample's mean seeing of 1.4", so if this plot were 
homogenized to the SDSS seeing most of these arcs would fall below our surface brightness 
limit. Seeing is thus the dominant limiting factor in our efficiency, and it is the relatively 
large PSF of the SDSS data that limits our detection to high surface brightness, quite long 
arcs with l/w > 10. On the other hand, the arcs that we can in fact flnd in the SDSS will 
therefore also be particularly interesting objects for follow-up observations. 
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Fig. 7. — The surface brightnesses of the arcs found in previous ground based searches. 
The surface brightness is reported in the AB system; given the simphfying assumption of a 
flat spectrum arc, we need not report the bandpass, though the data were in fact observed 
through a variety of filters {B,V, R, I). The vertical red line is the 50% efficiency limit of 
our visual search selection function. 
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6.1. Upper Limits on the Observable Arc Production 

We use our search resuhs to establish upper limits on the production of observable 
giant arcs by galaxy clusters. Assuming Poisson statistics, an observation of arcs for a 
given inspected cluster indicates that the 95% C.L. upper limit on the expected number 
of observable arcs per cluster is 3 (i.e., the probability of observing events for a Poisson 
distribution with mean 3 is 3°e~^/0! ~ 0.05 = 1 — 0.95). Then, given a null result after 
inspecting Ni clusters in a redshift bin i, the 95% C.L. upper limit on the expected number 
of observable arcs per cluster, Sg^^^, may be determined from the product of the individual 
Poisson probabilities using the relation 

[ (Sg^^j" / 0! = = 0.05 , (2) 

from which we find that ^ 

^95% = ~- (3) 



N, 



Sggy^ is shown in Fig. [H] for clusters with estimated masses M200 > 2 x lO^^M© and 
M200 > 4 X IO^^Mq. It should be noted that our upper limits correspond to arcs observable 
by the SDSS according to the efficiency for detection shown in Figs. [5] and [61 The lower 
panel in Fig. [8] corresponds to the more massive clusters (M200 > 4 x lO^^M©) in our sample, 
and from the plot it is clear that we only have the statistics for a significant limit in the 
range 0.1 < 2; < 0.3. 
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Fig. 8. — The 95% CL. upper limit for the probabihty for a cluster to produce a gravitational 
arc visible in SDSS imaging, assuming Poisson statistics in our sample of 825 clusters. The 
limit is shown as a function of redshift for M200 = 2 — 4 x lO^^Mo clusters in the top panel 
and for M200 > 4 x 10^"^ clusters in the bottom panel. 
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6.2. Comparison With Previous Works 



In order to understand the compatibility of our results with previous arc searches, we 
perform a direct comparison with these searches over similar redshifts. The cluster selection 
function is a considerable source of uncertainty which is hard to quantify for this comparison, 
ans so to reduce the impact of this effect we have selected on each sample clusters with 
similar mass estim ators. The sample of X-ray luminous galaxy clusters constructed by 
Smith et al.l (120051 ) is a good place to start. This sample has a total of 10 rich clusters with 
^^200 > 4 X 10^'^M (Tt at 0.17 < 2; < 0.26, where the mass is estimated from X-ray luminosities 
(j Allen et al.l 120031 ). Th e sample was imag ed with HST and the geometrical properties of the 
arcs were published by ISand et al.l (120051 ). In order to estimate the expected number of arcs 
observable in the SDSS, the parameters of the arcs are degraded from those observed in HST 
imaging (with a typical PSF of 0.15") to the image quality of the SDSS (with a typical PSF 
of 1.4"). This conversion is done assuming that the arc is well resolved in HST imaging, i.e., 
that the width is entirely due to the object, not the PSF. The width is then increased by 
adding in quadrature the average seeing for the SDSS, effectively reducing the l/w ratio. The 
increased width results in an increased area for the object, which translates into a reduction 
in surface brightness. Of the arc candidates observed with HST in these clusters, only 3 
remain with l/w > 10 after degradation to SDSS seeing. Of the surviving arcs only 1 has a 
surface brightness inside our detection limit. Thus the sample of ISmith et al.l (120051 ) gives a 
10% probability of having an arc observable by SDSS for clusters with M200 > 4 x lO^'^M© 
and 0.15 < z < 0.3. This result lies inside the 95% C.L. presented in Fig. [HI 



A similar comparison can be done with the iLuppino et al.l (119991 ) sample of 38 EMSS 
X-ray selected clusters with L^. > 2 x 10"^^ ergs/sec at 0.15 < 2; < 0.82. As before, the 
parameters for the arcs were transformed into expected p arameters for the SP SS using the 
seeing applicable to each arc observation carried out by iLuppino et al.l (119991 ) . After this 
procedure, no arc from this sample is expected to be observable in SDSS imaging. Again, 
this result is consistent with the upper limits derived from our analysis. 

The results o btained in this ^york ar e also consistent with a search for giant gravitational 
arcs in the RCS (iGladders et al.ll2003l ). The RCS cluster sample covers a redshift range 
0.35 < z < 0.95, and after visual inspection of approximately 900 clusters with M200 > 
4 X IO^^Mq, 3 strong lensing systems with giant arcs were found, all of them with z > 0.60. 
The result from the RCS indicates a smaller probability for a low redshift cluster to produce 
an arc, compared with a cluster of similar mass at higher redshifts. The bulk of the clusters 
used for the arc search in the present work have redshifts 0.1 < z < 0.3, lower than the arc 
clusters found by the RCS, and so our null result is also consistent with the RCS result that 
low-redshift clusters are less efficient for the production of arcs. 
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We are aware of 4 giant arc systenas reco vered or discovered in sin gle-pass SPSS imag ing;: 
one found by the RCS (iGladders et al.ll2003l . see below), one found by lSand et al.l (l2004j . RX 
J1 133), a new systern described in the Appendix, and another new system that is described 
in lAUam et al.l (120061 ) . Three of the arc systems are lensed by clusters with A''^ ~ 10 — 20, 
beneath the threshold of our current work, and these also have 0.35 < z < 0.45, beyond our 
catalog redshift limit though well within the reach of SDSS cluster finders. The fourth is 
RCS1419. 2+5326 at z = 0.64: our cluster finders can discover objects with these properties, 
but with Ng 5 (as observed by the SDSS; the real Ng is much higher) the catalogs are 
far from the 90% efficiency level of the z < 0.3 catalogs (the actual efficiency may be 
of order 1 — 5%). It is w orth noting that the two high surface brightness arcs found by 



Lynds &: PetrosianI (119891 ) are just south of the limits of the SDSS Stripe 82 region. Had 



both clusters been 1 degree farther north, they would likely have been found in our survey. 



7. Summary 

The potential for observing arcs in SDSS imaging has been demonstrated by the ob- 
servation of known arc systems in the mu ltiband imaging of the survey (see Fig. [1]) and by 



the discovery of new arcs (this paper and lAUam et al.ll2006l ). In this work, we presented a 



systematic search for gravitational arcs in a sample of SDSS clusters. A total 825 clusters 
with Ng > 20 (corresponding to M > 2 x IO^^Mq) have been visually inspected for arcs 
and have also been processed with an automated arc search algorithm. The efficiency of 
arc detection for both techniques has been quantified by using simulated arcs added to the 
inspected images. The results indicate a significantly higher efficiency for the eye scan and 
point to the need for developing a more efficient algorithm for the automated search. No 
gravitational arcs were seen in the sample of clusters analyzed. This result is consistent with 
previous surveys done over a similar redshift range (0.1 < z < 0.3). 
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A. Appendix: Hall's Arc 

We report here an arc discovered serendiptously by one of us (P.B.H) in 2004. It is 
widely known inside the SDSS Collaboration but has not been published previously. The 
purely serendipitous discovery was made during inspection of the images of spectra classified 
UNKNOWN (the initial reduction of the SDSS spectrum had bad spectrophotometry and 
could not be classified). There is a cluster clearly visible in the images: the BCG has a SDSS 
spectrum and is at z = 0.44, and it is coincident with NVSS J014656-092952. The SDSS 
spectrum shows two objects superimposed, the z = 0.44 galaxy and a star - the latter is 
what the SDSS pipeline reports. 

We ran the maxBCG code at the position of the BCG and measured Ng = 12, but this 
is a lower limit. The cluster is above the redshift where the Ng is complete - at this redshift 
two effects compromise the A^^ measurement: (1) 0.4L* is fainter than the useful limiting 
magnitude of the SDSS, and (2) the colors of objects near 0.4L* become noisy enough to 
scatter outside the color box used to determine cluster membership. A correction must be 
applied to Ng for these effects. As we are not attempting to use this arc in our analysis 
we defer the measurement of the corrected Ng to a later paper. We have tabulated our 
information about the cluster in Table [2l 

There are three arcs. The longest one was split by the SDSS photometry pipeline Photo 
into two objects, a and b. We report on the object parameters in Table [31 Three of the four 
objects have photometric redshifts of about 0.6-0.7; the fourth is about 3cr away (though we 
note that for faint objects r > 20 the SDSS photometric redshifts are less reliable because 
the spectroscopic training sets available are less extensive than at brighter magnitudes). We 
measured the parameters of the three arcs in the SDSS g-hand image. These are listed in 
Table HI It is likely that we would have found this arc in our survey (if it had been in our 
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sample catalog) because it is remarkably clean: there are no stars or cluster galaxies that are 
projected onto t he in iages of the arcs. Deeper imaging data for this system may be found in 



Hennawi et al.l (120061 ) . 
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Table 2. 


Hall's Arc: the Cluster 




Name 


RA 


Dec z 




SDSS+26.733-9.497+0.44 


01:46:56.00 


-09:29:52.4 0.44 


> 12^ 



Lower limit; see text. 
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Table 3. Hall's Arc: Positions, Magnitudes, and Photo-z's 



Arc 


RA 


Dec 


9 


u- g 


g-r 


r — i 


i — z 


^phot 


la 


26.73210 


-9.50009 


21.85 ±0.15 


1.3 ± 1.4 


0.57 ±0.21 


0.68 ±0.19 


0.09 ±0.51 


0.61 ±0.14 


lb 


26.73051 


-9.49967 


21.81 ±0.21 


0.39 ± 0.92 


0.75 ±0.27 


0.63 ± 0.22 


-1.8 ±2.6 


0.57 ±0.23 


2 


26.73063 


-9.49525 


21.66 ±0.10 


0.45 ± 0.46 


0.53 ±0.14 


0.73 ±0.12 


0.29 ±0.28 


0.67 ±0.08 


3 


26.73641 


-9.49673 


21.77 ±0.15 


0.66 ±0.81 


0.03 ± 0.27 


0.15 ±0.39 


-1.1 ±2.1 


0.15 ±0.15 



Note. — SDSS data. We retain the SDSS decimal degrees convention f or RA and Dec. The magnitudes and 
colors are SDSS model magnitudes and are corrected for extinction (see lAdelman-McCarthy et al.l l2006l . and 
references therein). 



'^The mean of the SDSS DR5 photoz and photoz2 (jAdelman-McCarthv et al.l 120061 ) photometric redshifts 
and errors. 
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Table 4. Hall's Arc: Arc Parameters 



Arc / l/w mean SB (mag/arcsec^) 



1 


17.2" 


11.2 


24.5 


2 


10.0" 


6.5 


24.6 


3 


3.8" 


2.5 


23.7 



Note. — SDSS (y'-band measurements. The 
5-band PSF is 1.53". 
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Name 


KA 




z 






bDbS+2d9.o+27.2+U.ll 


239.034 


27.1808 


A 1 AA 

0.109 


A A 

94 


O 1 o 

31.2 


bDbh+230.D+27.7+U.08 


O O A n 

230. D 


OT T1 /I /I 

27.7144 


A A O O 

0.083 


O O 

88 


OO 1 

28.1 


bUbb+213.D-00.3+0.13 


213.618 


A OO A1 

-0.3301 


A 1 OO 

0.133 


A 

69 


1 T T 

17.7 


CT^CC 1 1 ITT 1 A 1A 

bUbb+1 17.7+17.7+0. 19 


1 1 T T01 

117.721 


1 T d'^G O 

17.6768 


A 1 AO 

0.192 


68 


1 T 

16.7 


DUDk3+258.2+D4.0+0.08 


OCO OOT 


63.9924 


0.081 


67 


17.2 


CT^CC 1 OOT O 1 Ac: O 1 A AO 

bDbb+227.8+05.8+0.08 


OOT Tc: 

in. lb 


r TOOO 

5.7828 


A AO 1 

0.081 


65 


16.3 


C"r\CC 1 1 AT A A1 O 1 A OA 

oUlbb+197. 9-01.0+0.20 


1 AT O O/^ 

197.880 


1 OOOA 

-1.3329 


A OAO 

0.203 


65 


1 O 

15.3 


bDbb+2d4.9+34.4+0.25 


OO /I OAO 

234.893 


O /I /I O A 

34.4309 


A O /I A 

0.249 


O 

63 


1/11 

14.1 


bDbb+250. 1+40.7+0.25 


OCTA AAO 

250.098 


/IP TAOO 

40.7028 


A O /I O 


/^O 

62 


1 O T 

13.7 


bUbb+1 10.4+30.7+0. 15 


1 1 A O/^ 

110.30 


O TO O O 

30.7383 


A 1 /I O 

0.148 


O 

62 


1 /I 

14.4 


bDbb+227.6+33.5+0.12 


OOT o /I 

227.584 


O O 1 O 

33.486 


0.116 


A 

59 


1 O A 

13.4 


O T^O O 1 OCTAO 1 10 1 AOA 

bUbb+250. 8+13.4+0. 20 


OCA OAO 

250.843 


1 o O/^ o 

13.363 


A 1 AA 

0.199 


CT O 

58 


IOC 

12.5 


oJJoo-|-izD.o-|-4 / .i-|-U.io 


iZD.oZ 


A7 ^ AOA 


u.izy 


Oo 


1 Q 

iz.y 


SDSS+203.8+41.0+0.26 


203.818 


41.0132 


0.264 


55 


10.9 


SDSS+139.5+51. 7+0.24 


139.481 


51.7163 


0.236 


54 


10.7 


SDSS+255. 7+34.1+0.11 


255.09 


34.0611 


0.109 


54 


11.5 


SDSS+228.8+04.4+0.11 


228.823 


4.397 


0.107 


53 


11.1 


SDSS+137.3+11.0+0.18 


137.329 


10.9857 


0.179 


53 


10.7 


SDSS+184.4+03.6+0.08 


184.381 


3.6157 


0.089 


52 


10.9 


SDSS+216.5+37.8+0.16 


216.475 


37.7915 


0.164 


50 


9.7 



Table 5:: List of the most 20 most massive inspected 
clusters in our sample. The full sample is available elec- 
tronically. 
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